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ABSTRACT 

Despite technical advances, the future of chromatin 
mapping studies requires an ability to draw accurate 
comparisons between different chromatin states to 
enhance our understanding of genome biology. In 
this study, we used matched chromatin prepar- 
ations to enable specific and accurate comparisons 
of Saccharomyces cerevisiae chromatin structures 
in the presence and absence of the co-repressor 
protein Tup1. Analysis of wild-type and tupl A chro- 
matin data sets revealed unique organizational 
themes relating to the function of Tup1. Regulatory 
regions bound by Tup1 assumed a distinct chroma- 
tin architecture composed of a wide nucleosome- 
depleted region, low occupancy/poorly positioned 
promoter nucleosomes, a larger number and wider 
distribution of transcription factor-binding sites and 
downstream genes with enhanced transcription 
plasticity. Regions of Tup1 -dependent chromatin 
structure were defined for the first time across the 
entire yeast genome and are shown to strongly over- 
lap with activity of the chromatin remodeler Isw2. 
Additionally, Tup1 -dependent chromatin structures 
are shown to relate to distinct biological processes 
and transcriptional states of regulated genes, 
including Tup1 stabilization of Minus 1 and Minus 
2 promoter nucleosomes at actively repressed 
genes. Together these results help to enhance our 
mechanistic understanding of Tup1 regulation of 
chromatin structure and gene expression. 

INTRODUCTION 

The organization of eukaryotic DNA into chromatin is 
known to have far-reaching implications on the accessibil- 
ity and functionality of genetic information stored in the 



nucleotide sequence (1-3). Eukaryotic genomes are pack- 
aged through the wrapping of 147-bp segments of DNA 
around histone proteins arranged in octamers known as 
nucleosomes (4,5). Multiple isoforms of the histone pro- 
teins can combine to form a number of distinct nucleo- 
some octamers (6-8). Additionally, histones can undergo 
an assortment of covalent modifications to further diver- 
sify the chromatin landscape (9). 

The coiling of each 147-bp DNA segment around a nu- 
cleosome spans approximately 1.7 helical turns and estab- 
lishes important topological relationships between the 
nucleotide sequence and its local chromatin environment 
(5,10). Nucleosome positioning relative to underlying 
DNA sequences controls the access of regulatory factors 
to their genomic-binding locations and helps to coordinate 
expression programs in vivo (11-13). Consequently, the str- 
ength of each nucleosome's interaction with its associated 
DNA sequence can influence gene expression (14,15). 

In yeast and other eukaryotes, nucleosome positioning 
is controlled by a combination of inputs. Rigid poly(A) 
tracts of DNA serve as nucleosome-excluding sequences 
and help direct the formation of nucleosome-depleted 
regions (NDRs) both in vivo and in vitro (15-19). 
Chromatin architecture is also believed to be influenced 
by statistical packing principles whereby well-positioned 
nucleosomes help organize the regular positioning of nu- 
cleosomes in adjacent large stretches of DNA (20-22). In 
addition to sequence-based considerations and packing 
rules, a number of protein complexes can actively 
position nucleosomes to regulate chromatin structure and 
transcription (2,23-25). The Tupl-Ssn6 corepressor in 
Saccharomyces cerevisiae is one such complex, soliciting 
several chromatin-mediated mechanisms to regulate gene 
expression in vivo (26). 

In yeast, Tupl-Ssn6 is a global repressor of transcrip- 
tion, responsible for the repression of more than 1 80 genes 
involved in diverse signaling pathways (27,28). Tupl-Ssn6 
was one of the first corepressor complexes to be identified 
(29), and has since served as a model for similarly 
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structured corepressor proteins in other eukaryotes, inclu- 
ding the Drosophila Groucho (Gro) protein, and the mam- 
malian transducin (3-like/transducin P-like related (TBL/ 
TBLR) proteins and transducin-like enhancer of split 
(TLE) proteins (26,27,30). The Tupl-Ssn6 corepressor 
complex does not bind directly to DNA, but instead is 
recruited to promoters by sequence-specific DNA- 
binding proteins which coordinate the expression of specif- 
ic subsets of genes (27). Upon recruitment to promoters, 
Tupl-Ssn6 is known to repress downstream genes 
through several mechanisms, including interfering with 
the recruitment of transcription machinery, histone 
deacetylase (HDAC) recruitment, and the establishment 
of nucleosome positioning (26). Chromatin-dependent 
mechanisms of Tupl repression (nucleosome positioning 
and HDAC recruitment) are not mutually exclusive and 
redundant repression mechanisms are observed at many 
Tupl -regulated promoters (31). Moreover, unequal 
sensitivities of Tupl -regulated genes to the inactivation 
of different repression pathways in vivo suggests that re- 
sponses to Tupl repression may vary in a gene-specific 
fashion, whereby different groups of genes have evolved 
different strategies for utilizing Tupl-Ssn6 repression 
mechanisms. Genes that have developed a strong reliance 
on Tupl -dependent nucleosome positioning for repression 
in vivo have been well characterized by site-specific studies. 
Of the more than 100 genes regulated by Tupl only a few 
have been characterized in detail including RNR3 (32,33), 
FLOl (34), ANB1 (35), SUC2 (36,37) and other a-cell- 
specific genes (38^10). Tupl -dependent nucleosome pos- 
itioning at the promoter regions of these genes is believed 
to repress transcription by limiting the accessibility of 
promoter elements to trans-acting factors (26). 

In this study, we defined the genomic landscape for 
Tupl's regulation of chromatin structure, mapping chro- 
matin in both wild-type and tupl A cells using a combin- 
ation of matched micrococcal nuclease (MNase) digestions 
and high-throughput DNA sequencing methods. We 
identified Tupl -specific alterations in chromatin architec- 
ture at 96 genes across the genome and determined that 
Tupl stabilized the positioning of the —1 and —2 
promoter nucleosomes of these genes when Tupl is ac- 
tively repressing transcription. We determined that Tupl 
regulates chromatin organization at the majority of its 
targets by cooperation with the ATP-dependent chroma- 
tin remodeler Isw2. In addition, we also identified distinct 
chromatin architecture at Tupl -bound promoters consist- 
ing of a wide NDR with low occupancy/poorly positioned 
promoter nucleosomes, and a larger number and wider 
distribution of transcription factor-binding sites (TFBS). 
These distinctive Tupl -bound promoters regulate genes 
which have high transcription plasticity, suggesting a func- 
tional role for the Tupl-Ssn6 complex in regulating genes 
with varied expression. 

MATERIALS AND METHODS 

Yeast strains and growth conditions 

Wild-type (BY4741) and tupl A (BY4741/ tupl A ::KariMX) 
strains (Open Biosystems) were used for nucleosome 



mapping experiments. Tupl deletion was confirmed by 
PCR analysis of genomic DNA. Chromatin mapping ex- 
periments sampled from the genomes of yeast strains 
growing at 30°C in exponential phase (0.6-0.8 OD 600 ), 
as described previously (41). 

Analysis of Tupl ChlP-chip data 

To compare low-resolution (0.5-1 kb) ChlP-chip experi- 
ments of Hanlon et al. (42) to our high-resolution nucleo- 
some mapping data, ChlP-chip data were converted into a 
continuous binding profile. Both continuous binding (log 2 
ratio) and P-value profiles were generated by ChlPOTle's 
sliding window approach (43). The probability of Tupl 
binding to any gene's regulatory region was assigned as 
the minimal _P-value for the ChlPOTle windows overlapp- 
ing that region. In this approach, every gene is assigned a 
P- value for Tupl binding. Regulatory regions with a 
Bonferroni corrected P < 1 x 10~ 6 and mock uncorrected 
P> 0.001 in any growth condition were designated as 
strongly bound. Weakly bound regions were also selected 
as regions with a Boniforroni corrected P<lx 10~ 3 and 
mock uncorrected P> 0.001 (Table 1). 

Matched MNase digestions 

Preparation of chromatin DNA samples for sequencing 
was conducted as described previously (16,25,44), with 
several key revisions (Supplementary Figure SI A). 
Cross-linked chromatin-DNA complexes were extracted 
by mechanical disruption (bead beating; 4 x 1 min 
sessions, 2 min on ice between sessions). MNase titrations 
were carried out as described previously (25); however, an 
added level of standardization was achieved by quantify- 
ing the input of extracted chromatin loaded into each 
digest reaction; 1 mg total protein from whole-cell 
extracts per 200 ul digest reaction. Total protein concen- 
tration for whole cell extracts (chromatin extracts) was 
quantified using a Bradford assay (OD 595 ). Following 
MNase treatment and cross-link reversal, DNA was phe- 
nobchloroform precipitated, resuspended in 50 ul TE pH 
8.0, and RNase treated. Small aliquots (20%, ~10ul) 
from each digest titration were then analyzed by gel elec- 
trophoresis (Supplementary Figure SIB). Gel intensity 
measurements for each lane were calculated using standard 
densitometry software provided by Biorad (Quantity 
One™) and exported to Microsoft Excel™ for correl- 



ation analysis (Pearson). 


Correlation 


coefficients 


were 


Table 1. Tupl binding at yeast regulatory regions 




Regulatory region 


Abbreviation 


ChlPOTle. 


No. 






/■-value* 




Not bound 


NB 


>0.001 


5323 


Strongly bound unidirectional 


SBU 


<0.000001 


238 


Weakly bound unidirectional 


WBU 


<0.001 


189 


Strongly bound bidirectional 


SBB 


<0.000001 


160 


Weakly bound bidirectional 


WBB 


<0.001 


74 


Unknown 


Unk 


NA 


379 



*ChIPOTle /-values are corrected for multiple testing. Unknown is 
caused by bad IP data. 
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calculated for digest titration levels no longer showing 
a visible di-nucleosome band (i.e. optimally digested sam- 
ples) and spanning a region known to cover a size range of 
0-400 bp of DNA (Supplementary Figure SIC). A 100-bp 
standard ladder (Qiagen) was analyzed on the same gel to 
calculate a relative front (i.e. separation relative to 
standard fragment sizes) for the migrating DNA popula- 
tion. 'Matched' digests were selected as two optimally di- 
gested samples having a correlation coefficient (?->0.9). 
Correlation coefficients r > 0.9 were selected because 
these showed the most reproducible change in MNase pro- 
tection data on single site real-time quantitative PCR 
(RT-qPCR) analysis (See Supplementary Methods section). 
The remaining sample from each optimal digest was 
column purified (Zymo Research), bypassing the need to 
gel excise mono-nucleosome DNA. Avoidance of gel ex- 
traction to isolate chromatin DNA at this stage is a key 
step to assure sampling of similar chromatin populations. 
Matched DNA samples were prepared for sequencing 
using the standard Illumina protocol with additional care 
at the gel purification step to ensure that the same size 
range of DNA was selected. 

Nucleosome mapping 

Nucleosome DNA was sequenced by an Illumina Genome 
Analyzer IIx, as described previously (16,25,45). 
Sequencing reads were aligned to the March 2010 build 
of the 5. cerevisiae genome (www.yeastgenome.org) using 
the ELAND algorithm, allowing only two possible mis- 
matches. MNase protection was calculated for sequencing 
data and all comparison data sets by extending aligned 
reads by 120 bp, dividing by the average read count per 
base pair and converting this ratio into continuous space 
(log 2 ratio). Each data set was then standardized at single 
base pair resolution to have a mean occupancy of 0 and a 
standard deviation of 1 . Individual nucleosome calls were 
made using the template filtering algorithm with default 
parameters, the seven standard defined templates and a 
minimum and maximum allowable nucleosome width of 
80 and 220 bp, respectively (25). 

Nucleosome statistics 

Transcription start-site aligned data comparisons. Yeast 
transcription start sites (TSSs) were taken from David 
et al. (46). Significant differences for Tup 1 -bound and un- 
bound promoters between wild-type and tuplA average 
MNase protection plots were tested at 10-bp resolution 
for regulatory regions spanning TSS (—1000 to +200 bp) 
using a two-tailed ?-test. The resulting P-values were 
corrected for multiple comparisons by conservative 
Bonferroni correction. Tup 1 -bound promoters were 
sorted by downstream expression changes following 
Tupl deletion using data from Green et al. (28). Genes 
with decreased downstream gene expression following 
Tupl deletion were selected as being in the bottom 15% 
(15th percentile) of genes expressed in tupl A. Likewise, 
genes with increased downstream gene expression follow- 
ing Tupl deletion were selected as being in the top 15% 
(85th percentile) of genes expressed in tupl A. 



Template filtering comparisons. Nucleosomes called by 
template filtering were annotated to their closest promoter 
according to distance from TSS. The closest nucleosome 
upstream of each TSS, dyad within 500 bp, was annotated 
as —1, with the second closest as —2, and third as —3. 
Similar annotations were made for downstream +1, +2 
and +3 nucleosomes, respectively, for each promoter. 
Comparisons between called nucleosomes from unbound 
(Class NB) versus bound (Class SBU) promoters were 
made using 10000 permutation tests. Permutation 
testing randomly sampled the unbound promoter nucleo- 
some data set for the parameters under examination. A 
similar analysis was conducted when comparing bound/ 
wild-type versus bound/?M/?/zl promoter nucleosomes. 

RESULTS 

Matched MNase digests produce highly correlated 
chromatin samples 

To dissect the influence of Tupl on chromatin structure 
in vivo, we employed a technical approach which minim- 
izes the technical variation typically associated with 
MNase digestion steps utilized by genome-wide chromatin 
DNA preparations (see 'Materials and Methods' section; 
Supplementary Figure S1A) (25). Our approach allowed 
us to systematically select similarly digested chromatin 
preparations (i.e. matched MNase digestions) to control 
for differences in nuclease digestion extent and collect 
highly comparable chromatin populations. Matched 
digests were collected, similar to other chromatin-mapping 
protocols, but with the addition of several key steps (See 
'Materials and Methods' section) to help ensure greater 
consistency of DNA sampling. Overall nuclease digestion 
profiles were compared by analyzing chromatin DNA 
populations (i.e. DNA remaining following MNase treat- 
ment) separated by gel electrophoresis. Matched samples 
were selected as having a highly correlated (;• > 0.9) distri- 
bution of chromatin DNA population size and abundance, 
as identified by ethidium bromide DNA population staining 
intensities (Supplementary Figure SIB and SIC). Matched 
samples were then deeply sequenced in parallel using an 
Illumina Genome Analyzer IIx. 

Sequencing of matched chromatin samples generated an 
extensive sequencing data set for the two strains with 20.6 
and 23.7 million mapped sequencing reads for wild-type 
(BY4741) and tupl A genomes, respectively. Given that the 
yeast genome is packaged by an estimated 70 000 nucleo- 
somes (47), our data represents greater than 250 x average 
coverage for nucleosomal DNA. To our knowledge, 
the magnitude of this coverage provides the most exten- 
sive mononucleosome sequencing data set available for 
S. cerevisiae to date. Sequencing data was processed in 
three ways depending on the downstream analysis em- 
ployed. First, we examined the ends of the sequence 
reads on individual strands (Figure 1A). This representa- 
tion of the data is completely unprocessed and is an exact 
representation of the raw aligned sequencing output. 
Second, we extended each 36-bp sequence read to a total 
length of 120 bp. Although nucleosomes are 147 bp, we 
found that extrapolation to 120 bp provided better 
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Figure 1. Matched MNase digests identify well-characterized Tup 1 -dependent changes in chromatin structure at the FLOl promoter. Graphs 
illustrating sequencing data output, aligned to FLOl ATG for wild-type (Black) and tuplA (Red) experiments. (A) Raw sequencing data output 
(MNase Cleavage Density) displaying the number of sequence read ends which map uniquely to the forward (upward scale) and reverse (downward 
scale) DNA strands. (B) Mapped sequencing reads are determined by extending the uniquely mapped tags to 120 bp and standardizing the total tag 
count to 1 million per experiment. (C) MNase Protection is the log 2 ratio of the number of extended tags divided by the average read count per bp. 
Gray boxes: Regions of significant Tupl enrichment by ChIP (Bonferroni corrected P <1 x 1(T 6 ) (42). 



demarcation of nucleosome ends than did extrapolation to 
147 bp with no effect on occupancy comparisons. We then 
counted the number of times each unique position in 
the genome was included within a 120 bp extended se- 
quence and standardized between experiments to a total 
of 1 million sequence reads, to correct for variable se- 
quencing coverage (Figure IB). Lastly, extended read 
counts were converted to a 'MNase protection' score by 
dividing by the average tag count per base pair, and con- 
verting into continuous data-space by calculating a log 2 
ratio (Figure 1C), similar to previous analysis (16,45). We 
prefer to use the term 'MNase protection' to describe this 
data rather than 'nucleosome occupancy' because MNase 
protection may alternatively be conferred by non- 
nucleosomal protein complexes. MNase protection reflects 
the relative level of DNA protection from MNase diges- 
tion and at most locations represent the occupancy/ 
density of nucleosomes at each genomic location in rela- 
tion to the entire genome sequenced in that chromatin 
preparation (48). 

To test the homology between our two MNase protec- 
tion data sets, a Pearson correlation coefficient was calcul- 
ated comparing wild-type and tuplA chromatin maps 
across the entire genome. As expected, the correlations be- 
tween our matched samples was extremely high (r = 0.934). 
To compare this homology to other data sets, the same 
correlation analysis was repeated for other published 
chromatin data sets of wild-type yeast of the same back- 
ground genotype also grown in rich media. Across eight ex- 
periments from two different groups, only two experiments 



showed a correlation similar to ours (Supplementary 
Table SI). Interestingly, our matched samples, derived 
from different strains, were more similar than the majority 
of published replicates from the same strain. This result 
highlights how differences in technical preparation can 
introduce deviations into chromatin maps and how our 
matched preparation helps to control this influence. 

To better visualize local versus global homology be- 
tween our chromatin data sets, we calculated correlation 
coefficients and average changes in MNase protection for 
all 1-kb windows covering the yeast genome. As expected 
from our almost identical data sets, >90% of 1-kb win- 
dows had a correlation >0.9 and the majority of windows 
showed little change in average MNase protection (90% 
were within ±0.1) (Supplementary Figure S2). Overall, 
results from both approaches (correlation and average 
change) indicate that only a small fraction of the genome's 
chromatin changes in the absence of Tupl. 

As an alternative approach to assess the similarity of 
our samples, we utilized a computational algorithm for 
characterizing nucleosome DNA signals from high- 
throughput DNA sequencing methods, known as template 
filtering (25). We employed this approach with the seven 
available templates to determine locations and character- 
istics for all nucleosomes in both data sets, mapping 
57469 and 57 327 nucleosomes each for wild-type and 
tupl A chromatin, respectively (Supplementary Table S2). 
Importantly, analysis of nucleosomes called by template 
filtering confirmed that wild-type and tuplA chromatin 
samples were of similar size populations and had similar 
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template frequencies, as expected of a matched digest 
(Supplementary Figure S3). 

Importantly, our sequencing results also recapitulated 
results of published site-specific studies of nucleosome 
positioning at the PH05 promoter (Supplementary 
Figure S4) (13). Nucleosomes at the PH05 promoter 
have not been shown to be influenced by Tupl regulatory 
mechanisms and, as expected, homologous levels of 
MNase protection are seen at this locus (between wild- 
type and tupl A chromatin) compared to known Tupl- 
regulated regions including FLOl, SUC2 and RNR3, 
which showed significant dissimilarity that was specific 
to the loss of Tupl (See next section). 

Matched MNase digests identify Tupl -specific differences 
in chromatin structure 

Having successfully prepared highly correlated chromatin 
data sets, we hypothesized that the high degree of hom- 
ology between our samples would enhance our ability of 
detecting Tupl -specific differences in chromatin structure, 
when comparing our wild-type and tupl A data sets. To 
test this hypothesis, we explored the relationship between 
Tupl binding and chromatin structure differences between 
our data sets by integrating our sequencing data with a 
recently published ChlP-chip data set on Tupl binding 
(42). Specific genomic fragments were scored for Tupl 
binding with a P- value and/or Tupl occupancy (binding) 
score (see 'Materials and Methods' section). Depending 
on the analysis, 100-bp genomic segments or 1-kb regula- 
tory regions were scored. A gene's regulatory region was 
defined as extending 1-kb upstream of the open reading 
frame, and each such regulatory region was classified with 
respect to Tupl binding by P-value (Table 1). Each bound 
regulatory region was further classified as bidirectional or 
unidirectional depending on if it was upstream of two di- 
vergently transcribed genes (list of all regulatory regions 
available in Supplementary Table S3). For the majority of 
our analysis, we compared regulatory regions classified 
as not bound by Tupl (NB) to strongly bound unidirec- 
tional (SBU) and/or weakly bound unidirectional (WBU) 
(Table 1). Importantly, integration of genome-wide Tupl- 
binding data with our chromatin sequencing data sets 
revealed that our approach accurately reproduced the re- 
sults of many high-resolution single-site studies on the per- 
turbation of chromatin by Tupl. For example, at the 
well-studied Tupl -regulated locus, FLOl, peaks of Tupl 
enrichment aligned with close proximity to regions of 
chromatin dissimilarity between matched-digest wild- 
type and tupl A samples, consistent with the effects of a 
known direct-acting Tupl-Ssn6 repression mechanism at 
that location (Figure 1). Moreover, previously mapped 
differences in FLOl chromatin structure between wild- 
type and ssn6A strains were captured by our analysis with 
high fidelity and sensitivity, including changes in nucleo- 
some occupancy and positioning directly upstream from 
transcription start (—1000 to — 1 bp) and long-range dif- 
ferences (-5000 to -4000 bp) (34). In addition to FLOl, 
the well-documented role of Tupl in establishing nucleo- 
some positioning and occupancy at the SUC2 and RNR3 



promoters was also evident in our data (Supplementary 
Figure S5 and S6) (32,33). 

To examine our MNase protection data genome-wide in 
relation to Tupl binding, we plotted the MNase protec- 
tion profiles (i.e. nucleosome occupancy) for both data 
sets against each other for consecutive non-overlapping 
100-bp windows covering the entire yeast genome 
(Figure 2A). As expected from the high correlation coef- 
ficient for our data sets, the plot illustrated clearly that the 
majority of the genome's 100-bp segments maintained 
similar MNase protection values regardless of the presence 
of Tupl (R 2 = 0.86). Interestingly, however, regions with 
dissimilar MNase protection profiles between wild-type 
and tuplA data sets occurred in a directional fashion, 
whereby decreased nuclease protection was only seen for 
tupl A data relative to wild-type. The loss of MNase pro- 
tection (and underlying nucleosome occupancy) when 
Tupl is deleted is consistent with Tupl-Ssn6's role in 
stabilizing nucleosomes in vivo (26) (41). The specificity 
of this result can be discerned when plotting only regions 
of the genome containing the lowest Tupl -binding scores 
(Bottom 10%) and comparing to regions with the highest 
Tupl -binding scores (Top 10%) (Figure 2B and C, re- 
spectively). When the bottom 10% of Tupl -binding sites 
are plotted there is an absence of locations showing a loss 
in MNase protection when Tupl is deleted, and the simi- 
larity between chromatin data sets improves (R = 0.92). 
On the other hand, when examining the top 10% of 
Tupl -binding sites there is a large directional bulge and 
the chromatin data sets show increased dissimilarity 
(R 2 = 0.72). These data demonstrate that differences in 
MNase protection between our chromatin data sets, across 
the yeast genome, are specifically related to neighboring 
Tupl binding and occur in a directional fashion consistent 
with Tupl's role in stabilizing nucleosome occupancy. It is 
important to remember that the Tupl ChlP-chip data 
were obtained at low resolution (0.5-1 kb), and many of 
the sites with high Tupl -binding scores would be caused 
by proximity to true binding locations. Given this differ- 
ence in resolution, it is likely that the specificity of the plot 
in Figure 2C is underestimated. 

As an alternative approach to investigating the specifi- 
city of Tupl binding to chromatin structure differences 
between our data sets, regions of chromatin dissimilarity 
(1-kb windows or larger) were classified as either having a 
poor correlation between MNase protection data (r < 0.5) 
or a large average change in MNase protection across a 
given region of interest (tupA-wt < — 0.5). The average 
Tupl -binding score was calculated for these dissimilar re- 
gions and compared to a random sampling of similar-sized 
windows across the yeast genome (Figure 2D and E). 
Overall, regions of high dissimilarity in chromatin struc- 
ture were enriched for Tupl binding to a significantly 
greater extent than randomly sampled regions, indicating 
that high Tupl enrichment was a specific characteristic of 
dissimilar chromatin regions between our data sets, again 
consistent with Tupl's role in modifying chromatin struc- 
ture (26). Interestingly, the specificity of Tupl enrichment 
to chromatin dissimilarity between wild-type and tupl A 
data sets was eroded when unmatched wild-type chroma- 
tin preparations were instead compared to our tupl A data 
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Figure 2. Deletion of Tupl causes a loss in MNase protection at Tupl-bound sites. (A) Density plot displaying MNase protection profiles for 
wild-type versus tuplA data sets for all non-overlapping 100-bp segments spanning the entire yeast genome. (B) Density plot of 100-bp genome 
segments in 2 A which are not bound by Tupl [bottom 10% enrichment by ChIP (42)] (C) Density plot of 100-bp genome segments in 2 A that are 
bound by Tupl [top 10% enrichment by ChIP (42)]. (D) Average Tupl binding (log 2 ratio calculated from ChIP enrichment score, see 'Materials and 
Methods' section) seen at dissimilar regions, identified by correlation between wild-type and tupl A chromatin structures (Black bar) compared to 
Tupl enrichment seen at randomly sampled (n = 10000) genomic regions (Gray bar). Randomly sampled windows were composed of an identical 
window size composition to the dissimilar regions. Error bars represent standard error in average Tupl -binding scores across sampled windows. 
Dissimilar regions were defined as continuous windows (>lkb) with a correlation /-<0.5 between the data sets. (E) Same analysis as 2D for regions 
(>lkb) with a large average change in MNase protection between chromatin structures (tupl A-wt <— 0.5) compared to randomly sampled 
(it = 10000) genomic regions, again with an identical window size composition. 



(Supplementary Figure S7). The decrease in average Tupl 
enrichment seen at regions of dissimilarity between wild- 
type and tuplA chromatin structures when unmatched 
chromatin samples are compared highlights the utility of 
collecting highly correlated chromatin DNA populations. 

Tupl primarily influences chromatin structure at 
promoter regions 

Comparing wild-type and tupl A chromatin maps allowed 
us to dissect the genomic landscape of Tupl's influence on 
chromatin structure. In total, the dissimilar windows 
identified by analysis for Figure 2D and E cover <4% 
of the budding yeast genome sequence, consistent with 
Tupl's presumed role modifying select regions of chroma- 
tin, rather than more broadly perturbing the genome 
at-large. Of dissimilar regions, 78% of dissimilar regions 
were identified by both correlation and average change 
analysis, suggesting that the loss of Tupl regulation influ- 
ences both nucleosome occupancy and positioning in vivo. 



Localization of Tupl-Ssn6 is known to vary from gene 
to gene and its influence on chromatin structure has been 
documented at both intergenic regions and overlying 
coding sequences (26). Of the 238 unidirectional regula- 
tory regions bound strongly by Tupl (Class SBU) in our 
experiments, 92 (38.7%) were identified as having at least 
500 bp of dissimilarity (by correlation, average change, or 
both) between wild-type and tupl A with 31 (13%) coding 
regions identified as being dissimilar. In total, 14 of the 31 
dissimilar coding regions identified were also associated 
with a dissimilar regulatory region, indicating that 
Tupl's activity on chromatin structure is predominantly 
based at promoter regions, a feature seen throughout our 
analysis. 

Tupl-bound genes show a distinct open chromatin 
architecture at promoter regions 

To assess differences in promoter packaging within and 
between our chromatin data sets, genes were aligned by 
TSS and averaged MNase protection was plotted at 10-bp 
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resolution for a 1-kb window up and downstream of the 
TSS, similar to previous studies (16,45). As expected, 
average aligned MNase protection data were extremely 
similar between wild-type and tuplA strains for all yeast 
TSS overall (Figure 3A). Like typical yeast promoters, 
both wild-type and tuplA averaged MNase protection 
(nucleosome occupancy) plots displayed a NDR directly 
upstream of the TSS flanked by well-positioned nucleo- 
somes, especially over the coding region. The same char- 
acteristic pattern of MNase protection was also seen for 
all promoters not bound by Tupl (Figure 3B). However, 
when promoters bound by Tupl were examined alone, a 
markedly different pattern of chromatin structure was 



seen. As illustrated by Figure 3C, genes with promoters 
bound by Tupl (Class SBU) exhibited an increased NDR 
width upstream of the TSS, and a significantly different 
pattern of nucleosome occupancy and positioning sur- 
rounded that region compared to unbound (Class NB) 
promoters (Figure 3B). The unique chromatin structure 
profile associated with Tupl -occupied promoters was 
not a product of our data normalization and transform- 
ation procedures, since the same result was seen when our 
data was viewed in non-continuous count space 
(Supplementary Figure SF8). 

As an alternate analysis of promoter structures, we also 
analyzed our chromatin data by individually calling 



ALL Promoter Types 



1 1 1 

. A 


i i i i i i 

'A - 




\j ^— wild-type 


i i i 


--tupl A 





03 



en 
O 



c 0 
u 

O -1 

CL 

<U 
rt3 



CT1 

TO 
s— 

<u 
> 

< 



Unbound Promoters 



1 1 1 1 

- B 


1 1 1 1 1 
11 / V V \# \/\_, 




wild-type 


1 1 1 1 


tupl A 

i i i i i 



Tupl -Bound Promoters 



1 1 

. c 


1 1 




1 1 1 1 


* * 




1 


wild-type 






tupl A 


1 1 




p<.05 
i i i i 



wild-type 



Unbound Promoters 



100 bp 




wild-type 



Tup1-Bound Promoters 

**(P = 0.0001) . 



, "(p = 0.0029) , 
'-2 *-1 
f- 



*+1 



*+2 



*+3 



H \ 



TSS^^ 



H h 



o (p = 0.0001) a (p < 0.0001) 

tup1A_ 2 



a (p = 0.0032) 



a (p = 0.0039) 



cf (p = 0.0103)' 
+1 +2 +3 




A (p< 0.0001) A (p< 0.0001 



Nucleosome Occupancy 



-1000 0 

Distance (bp) from TSS 



1000 



Figure 3. Tupl -regulated promoters have a more open chromatin structure compared to the rest of the yeast genome. (A-C) Graphs of average 
MNase protection profiles aligned by TSS and sorted into specific promoter classes (Table 1) for wild-type (Black) and tupl A (Red) data sets. Only 
genes with mapped TSS are shown (46). (A) Average profiles for all promoter types listed in Table 1, excluding those of unknown type (n = 4741) (B) 
Average profiles for unbound promoters (Class NB; n = 4298). (C) Average profiles for Tupl unidirectional-bound promoters (Class SBU; n = 162). 
Gray boxes: Regions of significant change in average MNase protection following Tupl deletion, as determined by two-tailed ^-testing (Bonferroni 
corrected / > <0.05). (D) Figure drawn to scale, depicting promoter nucleosomes for wild-type (Black ovals) and tupl A (Red ovals) chromatin data 
sets. Promoter nucleosome statistics were calculated for nucleosomes called by the template filtering algorithm (25) and annotated to individual 
promoters as described in 'Materials and Methods' section. Relative nucleosome occupancy (density at each locus) is represented by shading inside 
nucleosomes (White = low, Black = high). Asterisks (*) denote significant decrease in nucleosome occupancy compared to unbound promoters. 
Double asterisks denotes significant increase in distance (bp DNA) compared to unbound promoters. Sigma denotes significant increase in the 
variability of nucleosome positioning (i.e. the standard deviation for called nucleosomes dyad locations) compared to unbound promoters. Open 
triangle denotes a significant loss in nucleosomes called by template filtering compared to wild-type Tupl-bound promoters (Class SBU). 
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promoter nucleosomes using template filtering (25). This 
analysis revealed that the —2, —1, +1, +2 and +3 promoter 
nucleosomes all showed significantly lower nucleosome 
occupancy values and significantly more deviance in pos- 
itioning at Tup 1 -bound promoters, compared to their 
unbound counterparts (Figure 3D). As suggested by the 
broader NDR, the internucleosomal distance between the 
— 1 and +1 nucleosomes was significantly larger (40-bp 
increase; P = 0.0001). Additionally, a significantly incre- 
ased inter-nucleosomal distance between the —1 and —2 
nucleosomes (29-bp increase, P = 0.0029) was also evident 
compared to unbound promoters. Distances between +1, 
+2 and +3 nucleosomes remained unchanged. 

In addition to examining our MNase protection data in 
relation to Tupl binding, we also preformed unsupervised 
clustering to categorize promoters across the genome in 
relation to our two data sets (Supplementary Figure S9). 
A'-medoids clustering with values of K > 3 consistently 
produced a discrete cluster with a wide NDR region, simi- 
lar to that seen for Tupl -bound promoters in Figure 3C. 
Hypergeometric distribution analysis for this wide-NDR 
cluster consistently revealed an enrichment for Tupl- 
bound promoters and also enriched for promoters with 
a GO Slim annotation for genes involved in translation. 
These results suggest a possible link between the wide 
NDR structure we are seeing, Tupl binding, and specific 
biological pathways. 

Tupl-bound promoters are accessed by more transcription 
factors which bind across a wider region 

Having characterized a distinct chromatin structure asso- 
ciated with Tupl-bound promoters, we next sought to 
identify additional mechanisms that could best explain 
the observed pattern of MNase protection (and underlying 
nucleosome positioning/ occupancy) at these loci. 
Previous work has shown that both TFBS locations and 
active (occupied) binding sites in vivo are preferentially 
enriched at the characteristic NDR region upstream 
from each promoter's TSS, consistent with the model 
that nucleosome occupancy has a strong influence on tran- 
scription factor (TF) occupancy (49-53). Given these 
characteristic distributions and their relationship to chro- 
matin structure at promoter regions, we tested whether 
differences in TFBS locations underscored the distinct 
chromatin structure profiles associated with Tupl-bound 
and unbound promoters. Using the locations of highly 
conserved TFBSs identified by Maclsacc et al. (51), we 
plotted the distribution of conserved TFBS locations in 
relation to downstream-coding regions (ATG) and over- 
lying chromatin data (Figure 4A and B). As illustrated by 
Figure 4A and B, promoters not bound by Tupl showed a 
characteristic distribution of TFBS locations, with loca- 
tions correlating to overlying occupancy profiles and sig- 
nificantly enriching at the characteristic NDR region 
(—10- to — 200-bp upstream of ATG) upstream of 
coding regions, similar to previous analysis (49). In 
contrast, Tupl-bound promoters showed a wider distribu- 
tion of conserved TFBS locations, lacking the peak of 
enrichment of TFBS locations at the characteristic NDR 
region (—100 to —200 bp). Conserved TFBS locations at 



Tupl promoters were instead more spread out in their 
distribution across the promoter region, consistent with 
the distinct wide-NDR architecture of these genes. 
Importantly, the unique distribution of TFBS locations 
was also seen for conserved TFBSs active in the condition 
tested (bound in YPD), was observed when compared to 
wide unbound promoters, and was not dependent on pro- 
moter type (single versus double) or Tupl -binding class 
(Data not shown). 

Knowing that sequence-specific DNA-binding TFs 
must compete with and/or interact with nucleosomes to 
bind DNA and regulate transcription, we hypothesized 
that an increased number of TFs and/or TFBSs could 
be present at Tupl-bound promoters, either causing or 
resulting from the unique chromatin architecture and 
underlying TFBS distributions we identified. To test this 
hypothesis, we calculated the average number of con- 
served TFBSs and TFs present for each promoter class 
in Table 1. As illustrated by Figure 4C and D, Tupl- 
bound promoters showed a significantly larger number 
of both conserved TFBSs and bound TFs compared to 
unbound promoters, on average. Importantly, significant 
differences were seen regardless of Tupl -binding class and 
an even larger difference was evident for both TFBS 
number and number of bound TFs when we compared 
Tupl-bound single promoters versus unbound single pro- 
moters (data not shown). 

Tupl-regulated genes are more plastic in downstream 
expression 

Having uncovered a distinct structure for Tupl-bound 
genes, we hypothesized that the wider NDR with the 
higher number and spreading of conserved TFBSs would 
allow Tupl-regulated genes to respond more readily to 
various environmental signals. Accordingly, we proposed 
that Tupl binding would be associated with the capacity 
to modulate gene expression upon changing environmen- 
tal conditions, i.e. transcriptional plasticity (54). To 
examine transcriptional plasticity we used an extensive 
microarray data set for yeast grown in 13 diverse environ- 
mental transitions (55) and calculated transcriptional plas- 
ticity as the standard deviation for a gene's expression 
across all growth conditions. As illustrated in Figure 5A, 
our analysis found that Tupl binding was significantly 
correlated (P = 4.0 x 10~ 55 ) to transcriptional plasticity. 
Importantly, we also tested the specificity of the relation- 
ship between Tupl binding and downstream expression 
plasticity by conducting the same analysis for all regula- 
tory factors bound in rich media. Interestingly, we found 
that the association of Tupl -binding with downstream 
plasticity is stronger than that seen for all other factors 
tested, implying that plastic gene expression is a specific 
property of Tupl-regulated genes (Figure 5B). Overall, the 
strong association of Tupl binding with downstream ex- 
pression plasticity suggests that the distinctive packaging 
of Tupl-bound promoters may reflect a shared ability to 
broadly alter gene expression programs via Tupl regula- 
tion to coordinate appropriate responses to environmental 
signals. 
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Figure 4. Tup 1 -bound promoters show a unique distribution of TFBSs and are targeted by more transcription factors than unbound promoter 
regions. (A) Comparison of average MNase protection profiles (i.e. nucleosome occupancy) and (B) the underlying distribution of conserved TFBS 
locations at unbound (Class NB) and Tupl-bound (Class SBU) promoters, aligned by ATG. Data was plotted relative to translation start rather than 
transcription start due to the larger number of available TFBS locations and chromatin data using this approach. (C) Comparison of the average 
number of conserved TFBSs present at unbound promoters (Class NB) compared to Tupl-bound (Class SBU) promoters. (D) Comparison of the 
average number of TFs bound at unbound promoters (Class NB) compared to Tupl-bound (Class SBU) promoters. Double asterisks denotes a 
statistically significant increase (P < 0.0001) for (C) and (D). Error bars represent standard error in average signal detection. 



Tupl stabilizes the —1 and —2 promoter nucleosomes 

Having identified several structural features unique to 
Tupl-bound promoters, we next sought to examine how 
the deletion of Tupl influenced this chromatin structure, 
hoping to cast light onto the mechanism by which Tupl 
regulates chromatin. Analysis of our TSS plot of Tupl- 
bound promoters (Figure 3C) revealed that deletion of 
Tupl caused a decrease in MNase protection across the 
entire intergenic region with significant decreases at —1000 
to — 450-bp upstream of the TSS and directly upstream of 
the TSS; —70 to —10 bp (gray bars; Bonferroni corrected 
P < 0.05). Significant decreases in intergenic nuclease pro- 
tection levels overlied the positions of the —2 and —1 
promoter nucleosomes, upstream of the wider NDR, 
and suggest that Tupl -dependent nucleosome occupancy 
at these locations may play a key role in its regulatory 
mechanism. Additionally, when looking at individual nu- 
cleosomes called by template filtering, there was also a 



significant loss of called —1 and —2 promoter nucleosomes 
in the tupl A strain (Figure 3D). 

Significant changes in the average MNase protection 
profile in Figure 3C may appear slight because this plot 
encompasses all Tupl-bound promoters, including loci 
where Tupl is not regulating through chromatin. 
Therefore, to test the reproducibility of changes in 
MNase protection driving the significant differences 
identified in Figure 3C, changes in MNase protection 
and underlying histone H3 and H4 occupancy were 
assayed using RT-qPCR on three additional biological 
replicates at fifteen selected sites. As illustrated by 
Supplementary Figures SF10-SF12, strongly dissimilar 
regions between our data sets (r<0.5; YMR319C, 
YLR295C, YBR157C) showed large reproducible 
changes in MNase protection on RT-qPCR analysis and 
weakly dissimilar regions (0.5 < r < 0.7; YIL056W) 
showed smaller changes also with high consistency. 
Importantly, changes in MNase protection were not seen 
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Figure 5. Tupl binding at intergenic regions is specifically associated with transcription plasticity of downstream genes. (A) All genes for which there 
was data (n = 5802) have been sorted by Tupl occupancy upstream (42) and the moving average (50 genes) for Tupl occupancy versus expression 
plasticity is plotted. Expression plasticity was determined as the standard deviation across all 93 gene expression microarray experiments for various 
stress and metabolic responses (55). (B) Histogram displaying the correlation coefficients of genome-wide ChlP-chip data sets (204 proteins) versus 
expression plasticity (50,55). Tupl (indicated by the arrow) has the highest correlation to expression plasticity (r = 0.204; f = 4x 10~ 55 ). 



at the well-characterized PH05 locus, an expected result 
since this locus is not regulated by Tupl (Supplementary 
Figure S10). Additionally, significant changes in nucleo- 
some occupancy were also shown to underlie changes in 
nuclease protection profiles at all dissimilar regions tested 
thus demonstrating that at the majority of sites changes in 
MNase protection are caused by changes in nucleosome 
occupancy (Supplementary Figures S 1 O S 1 2) . Given these 
results, we were confident that strongly dissimilarity 
regions defined by an r < 0.5 were specific and reprodu- 
cible, therefore, the remainder of our analysis focused on 
strongly dissimilar regions between our two data sets. 

As an alternative approach to investigate Tupl stabil- 
ization of promoter nucleosomes in relation to transcrip- 
tional regulation, we again performed an analysis of 
TSS-aligned data for Tupl -bound promoters, but add- 
itionally sorted the chromatin data into three groupings 
according to downstream expression in tupl A strains (28) 
(Figure 6A-C). Genes repressed by Tupl and subsequent- 
ly expressed (de-repressed) upon its deletion, again ex- 
hibited a similar decreases in MNase protection across 
the regulatory region, significantly at the location of the 
—2 nucleosome (Figure 6B, gray bars; Bonferroni cor- 
rected P<0.05). Interestingly, an opposing result was 
obtained at genes where Tupl may be playing a role in 
activation of expression (i.e. repressed in tupl A). In place 
of losing nuclease protection at activated promoters of 
these regions, we instead saw a slight increase in protec- 
tion over the NDR region (Figure 6C), indicating that 
Tupl's influence on the chromatin at these regions is dif- 
ferent from that at genes it actively represses. 
Additionally, when examining genes with expression 
profiles unchanged by Tupl deletion in the condition 
being examined (rich medium; Figure 6A), no change in 
chromatin structure was seen, suggesting that the differ- 
ences we see in chromatin structure for Tupl -bound genes 
are due to its role in repressing transcription. 



In addition to sorting MNase protection data for 
Tupl -bound promoters by downstream expression, we 
clustered change in MNase protection profiles (tupl A — 
wt) for all Tupl -bound promoters for a 1-kb window 
surrounding TSS (Figure 6D). Clustering of bound pro- 
moters yielded only two significantly distinct patterns of 
change in MNase protection, regardless of the number of 
clusters generated. Cluster 1 had a broad increase in 
MNase protection across Tupl -bound promoters and 
consistently enriched for genes with associated GO slim 
process of translation, the same functional term associated 
with Tupl -bound promoters in our initial clustering 
analysis (Supplementary Figure S9). The broad increase 
in MNase protection seen at these promoters also parallels 
the increase seen for genes believed to be activated by 
Tupl (Figure 6C). Cluster 2 showed a broad decrease in 
MNase protection across promoter regions and signifi- 
cantly enriched for downstream genes de-repressed in a 
tupl A strain, similar to the result in Figure 6B. The 
same changes in structural patterns for Tupl -bound 
genes were also seen when we clustered change in 
MNase protection data for all available TSSs (including 
unbound promoters). As illustrated in Supplementary 
Figure S13, broad increases or decreases in MNase pro- 
tection are again seen for wide-NDR promoters that 
enrich for Tupl binding (Clusters 1 and 2). Similar to 
results in Figure 6D-F, Tupl -associated promoters with 
a broad increase in MNase protection following deletion 
of Tupl again enrich for genes associated with translation 
and ribosome biogenesis and assembly. In contrast, 
Tupl -associated promoters with a broad decrease in nu- 
cleosome occupancy following deletion of Tupl enrich for 
genes involved in response to chemical stimuli and 
unknown biological processes. Together, results from 
Figure 6, Supplementary Figures S9 and S13 illustrate 
how Tupl can have different effects on chromatin 
structure, depending on both the associated biological 
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Figure 6. Tup 1 -dependent chromatin structure changes depending on the co-repressor's role in transcriptional regulation. (A-C) Graphs of the 
average MNase protection profiles for Tupl-bound promoters (Class SBU), sorted according to downstream expression changes following Tupl 
deletion (see 'Materials and Methods' section) (28). Both wild-type (Black) and tupl A (Red) MNase protection data are plotted. Gray boxes: 
Regions of significant change in average MNase protection following Tupl deletion, as determined by two-tailed t- testing (Bonferroni corrected 
P < 0.05). Only a few regions of repressed promoters demonstrated a significant loss of MNase protection because of the limited number of 
promoters in each category and the use of multiple testing corrections. For this analysis we tested a 1-kb region upstream of TSS (—1000 to 
Obp) at 10-bp resolution. DS = Downstream. (D) Heat map of relative changes in MNase protection data (tupl A-wt), aligned to TSS, for all 
Tupl-bound promoters (combining Class SBB, SBU, WBB and WBU). All four Tupl-bound promoter classes were included in this analysis to 
increase the number of promoters with available TSS for clustering. Change in protection data was clustered using unsupervised A--mediod clustering 
(k = 2) for data spanning 1-kb window up and downstream of TSS. (E) Table of statistically significant (Bonferroni corrected P< 0.0001) enrich- 
ments of GO Slim terms associated with genes in each defined cluster. Significant enrichments were identified by hypergeometric testing. (F) 
Graphical illustration of the average MNase protection data for each cluster. (Black) wild-type and (Red) tupl A data sets. 



processes and downstream expression events involved at 
the genes it regulates. 

Tupl stabilizes nucleosomes at stress response TFBSs 

Chromatin-dependent mechanisms of Tupl repression are 
believed to repress transcription by limiting the accessi- 
bility of promoter elements to trans-acting factors (26). 
To explore the relationship between Tupl regulation, 
chromatin structure and rrans-activating factors, we ex- 
amined the chromatin structure surrounding all annotated 
TFBSs (51). The average change in MNase protection 
(tupl A-wt) for each TF was determined, and change- 
in-protection data (tuplA-wt) were aligned by TFBS 
center. To identify relationships among the TFs, TFBS 
change-in-protection data were sorted by the average 



change-in-protection for the middle 100-bp centering on 
the TFBS (Figure 7A). Three distinct patterns of chroma- 
tin structure emerged (Figure 7B-D) whereby the average 
change in MNase protection profiles increased, decreased 
or did not change following Tupl deletion (tupl A-wt). 
Overall, the majority of chromatin at TFBSs was not 
altered by Tupl deletion, consistent with the fact that 
Tupl regulates just a small fraction of the genome. 

TFBSs which are centered on locations where MNase 
protection was decreased, suggesting nucleosome loss in a 
tupl A, could represent TFs whose binding is blocked 
directly by Tupl -stabilized nucleosomes or blocked by 
Tupl itself. Indeed, several of the TFs (Migl, Skol, 
Nrgl, Phdl, Sutl, Aftl, Cin5, Skn7, Yap6, Rcsl, Roxl) 
that showed a dramatic loss of MNase protection in the 
tuplA strain were known co-factors of the Tupl-Ssn6 
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Figure 7. Tup 1 -dependent MNase protection overlies stress response factor binding sites, including many Tupl recruiters. (A) Heat map of average 
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factors (MIG1 — Green, GAT1 — Black, FLH1 — Red) representing the upper, middle, and lower percentiles of sorted data displayed in 7A. 



complex. These TFs have been shown to physically interact 
with Tupl or Ssn6 and are likely DNA-binding recruiting 
proteins for the Tupl-Ssn6 complex (42). Identification of 
these factors further confirms the specificity of our assay 
at detecting Tupl -specific changes in chromatin structure. 



Additionally, other TFs with TFBSs present at MNase 
protected regions dependent on Tupl may represent still 
uncharacterized Tupl-Ssn6 cofactors or TFs whose 
binding may be regulated by Tupl. The general stress 
TFs Msn2 and Msn4 fall in this category and showed 
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decreased nuclease protection (i.e. nucleosome loss) in 
tuplA. Interestingly, Msn2/Msn4 have been shown to 
regulate many of the genes required for the environmental 
stress response (56-58), and our findings suggest that their 
binding may by regulated by Tup 1 -dependent nucleosome 
occupancy/positioning. 

Surprisingly, there were a few TFBSs located at regions 
showing an increase in MNase protection following Tupl 
deletion, including Rapl, Fhll, Sfpl and Abfl. All four of 
these TFs target protein biosynthesis genes, some of the 
most highly expressed genes in rapid growing yeast. Since 
Tupl has been shown by ChlP-chip analysis to bind at 
many protein biosynthesis genes, this change in nuclease 
protection could be directly related to Tupl regulation 
(42). To test if this increase in MNase protection was a 
direct or indirect effect of Tupl regulation, we examined 
all genes annotated with the protein biosynthesis GO-slim 
term and compared between Tupl -bound and unbound 
regulatory regions (Supplementary Figure SI 4). Protein 
biosynthesis genes showed an increase of similar magni- 
tude in MNase protection in the NDR upstream of the 
TSS, regardless of Tupl -binding class, suggesting that the 
change in MNase protection at these sites was not directly 
caused by Tupl. While this increase in MNase protection 
appears indirect, the overall architecture of Tupl -bound 
protein biosynthesis genes was significantly different from 
that of unbound protein biosynthesis genes. As we demon- 
strated earlier for all Tupl -bound genes, the NDR for 
Tupl -bound protein biosynthesis genes was significantly 
wider than for unbound protein biosynthesis genes and 
promoter nucleosomes showed variable positioning. This 
result suggests that there is a unique subset of protein 
biosysnthesis genes which is regulated by Tupl. 

Isw2 remodels the majority of Tupl remodeled chromatin 

Having identified distinct effects of Tupl on chromatin 
structure that relate to the associated biological processes 
and expression state of regulated genes, we next asked if 
Tupl's effect on chromatin structure was also associated 
with the activity of other regulatory factors. Tupl-Ssn6 
has been shown to require the ATP-dependent chromatin 
remodeling enzyme Isw2 to position nucleosomes for gene 
repression at the well-characterized RNR3 locus, both 
upstream of and across the RNR3-coding region (32,33). 
Given the well-established collaboration between Isw2 
and Tupl nucleosome-positioning activities in vivo at 
RNR3 and the highly similar effects that Tupl and Isw2 
yield on upstream nucleosome occupancy (59-61), we 
hypothesized that the association of Isw2-remodeling 
with Tupl -remodeling events extended beyond individual 
loci and represented a more frequent collaboration of re- 
modeling activities across the genome. To test this hypoth- 
esis, association testing was performed comparing Isw2- 
remodeled genes defined by Whitehouse et al. and Tupl- 
remodeled promoters identified by our analysis. Only 5' 
Isw2 remodeling events were considered for this analysis, 
since the regions of significant Tupl remodeling, we 
identified were at the 5'-end of genes (i.e. within 0- to 
— 1000-bp upstream of TSS, Figure 3C). As illustrated in 
Figure 8 A, a greater proportion of Tupl -remodeled sites 



were also remodeled by Isw2 (52%), when compared to 
non-Tup 1 -remodeled genes (21%). Importantly, with a 
relaxed stringency for defining Tupl remodeling (r<0.7) 
and inclusion of all Tupl -bound regulatory regions (SBU, 
WBU, SBB and WBB), there were a total of 96 Tupl 
bound and remodeled genes of which 48 or 50% were 
also Isw2-remodeled. A significant association of Tupl 
remodeling with Isw2 remodeling was also seen at un- 
bound (Class NB) promoters and likely represent direct 
effects whose Tupl -binding state was not identified by 
previous low-resolution ChlP-chip experiments (42). 
Overall, the significant association of Tupl remodeling 
with Isw2 remodeling indicates that collaborative actives 
of Tupl and Isw2 seen at the RNR3 promoter extend 
beyond individual loci and represent a genome-wide regu- 
latory collaboration. 

Hypoacetylation of promoter histones is associated with 
active Tupl repression 

Decreased histone acetylation has been shown to track 
with Tupl-Ssn6 localization (62,63) and it is believed that 
Tupl may repress gene expression by creating a self- 
reinforcing state through its interactions with hypoacety- 
lated histones and histone deacetylase complexes (HDACs) 
in vivo (26,63-66). We hypothesized that Tupl remodeling 
could be associated with a distinct histone modification 
profile at remodeled promoters, reflecting potential collab- 
orative activities of Tupl chromatin regulation with 
HDAC activity to regulate transcription. To test this hy- 
pothesis, we aligned average chromatin modification data 
from Pokholok et al. (67) for a 1-kb window surrounding 
TSS and compared modification profiles by Tupl -binding 
class. Interestingly, no significant differences in overall 
chromatin modification profiles were evident when com- 
paring Tupl -bound promoters to unbound promoters 
(Figure 8B). Instead, only when Tupl-bound promoters 
were separated into Tupl repressed genes versus non- 
repressed was hypoacetylation shown to track with Tupl 
targets which are actively being repressed (Figure 8B). 
Additionally, a similar modification profile was also seen 
for Tupl-bound promoters that were either remodeled by 
Tupl or remodeled by Isw2 (Supplementary Figure SI 5). 
This result suggests that Tupl's interactions with 
hypoacetylated histones and HDACs in vivo is a feature 
of actively repressed Tupl-bound loci and not Tupl local- 
ization alone. Additionally, that the distinct chromatin 
modification profiles associated with gene repression mir- 
rored those of Tupl- and Isw2-remodeled loci, further 
suggests that collaboration of HDAC complexes and 
other chromatin modifiers is a general feature of active 
Tupl repression through chromatin. 

DISCUSSION 

Matched MNase digests reduce technical variation 
and reproducibly identify Tupl-specific changes in 
chromatin structure 

Genome-wide chromatin DNA isolation protocols typic- 
ally utilize a combination of nuclease digestion and gel 
electrophoresis to isolate mono-nucleosomal DNA 
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Figure 8. Chromatin remodeling by Tupl is strongly associated with Isw2 remodeling activity. (A) Pie graphs illustrating the association of 5' 
Isw2-remodeling events with Tupl -remodeling events (remodeled denned as r<0.5 for 1-kb window directly upstream from TSS). Graphs comparing 
Tupl-bound (Class SBU) promoters and unbound (Class NB) promoters. Association data was sorted by promoter class and Tupl -remodeling state. 
(Red box): Tupl -remodeled promoters. (Black box): Non Tupl -remodeled promoters. (Yellow): Percentage of promoters remodeled by Isw2. (Blue): 
Percentage of promoters not remodeled by Isw2. The 8% of unbound (NB) genes remodeled by Isw2 was expected and is consistent with the rate of 
genome wide 5' Isw2 remodeling identified by Whitehouse and colleagues (60). (B) Average histone modification profiles from Pokholok et al. (68) 
for 1-kb window up and downstream of TSS. Promoters are sorted as unbound by Tupl (Class NB) or Bound (combining Class SBB, SBU, WBB 
and WBU). All four Tupl-bound promoter classes were included in this analysis to increase the number of promoters with available TSS and 
modifications data. Additionally, Tupl-bound promoters are sorted according to downstream expression in tuplA strain (see 'Materials and 
Methods' section). 



fragments (44). However, since MNase digestion does not 
release all nucleosomes at the same rate (68), the resulting 
mononucleosomal DNA population has a variable com- 
position which depends largely on the extent of nuclease 
digestion (25). This variation presents a challenge of dis- 
tinguishing genuine heterogeneity in nucleosome positions 
and occupancy/density between two chromatin popula- 
tions from those created as an artifact of technical prep- 
aration, a dilemma increasingly encountered by chromatin 
biologists (69) and one made significantly more difficult 
when attempting comparisons between data sets from dif- 
ferent laboratory groups, strains or organisms. Indeed, the 
brief comparison of genome-wide data sets in 
Supplementary Table SI highlights the influence that 
subtle technical variations alone can exert between indi- 
vidual biological replicates from the same strain. While 
comparison of our tupl A data to samples from the 
Weiner et al. digestion series made in Supplementary 
Figure S7 represent intentional exaggerations of the influ- 
ence that technical variation can yield on signal specificity, 
the low Tupl enrichment (i.e. poor signal specificity) seen 
for all Weiner samples indicates that genome-wide data set 
correlations up to, and possibly exceeding, r = 0.77 are 
insufficient to achieve specific comparisons of chromatin 
structures. 



Ultimately, the future of chromatin mapping studies 
requires an ability to make sensitive and specific compari- 
sons between different chromatin states, be that strains, 
cell-types or environmental conditions, etc. to enhance our 
understanding of genome biology. Our matched MNase 
digestion approach was designed to accomplish this goal 
by reducing the technical variation introduced by the 
MNase digestion extent and gel extraction steps common 
in genome-wide chromatin mapping protocols. The pro- 
cedure was successful in generating nucleosome maps that 
were essentially identical for the wild-type and tuplA 
strains sequenced, except for Tupl -specific differences in 
chromatin structure (Figure 2). More importantly, the 
specific differences in MNase protection our analysis iden- 
tified were highly reproducible and related to changes in 
underlying nucleosome occupancy (Supplementary 
Figures S10-S12). Together, this data underscores the 
utility of matching the extent of MNase digestion as an 
approach to reduce technical variance in chromatin prep- 
arations and to enable accurate and reproducible down- 
stream comparisons of chromatin structures. 

Tupl regulates chromatin at stress response genes 

Comparison of TSS-aligned promoter MNase protection 
data indicates that Tupl's influence on chromatin 
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structure is primarily localized at upstream regulatory 
regions, specifically those associated with the Minus 1 
and Minus 2 promoter nucleosomes (Figure 3C and D). 
Analysis of gene expression data also suggests that 
Tupl promotes increased MNase protection at upstream 
locations to facilitate the repression of gene expression 
(Figure 6B). Additionally, analysis of clustering data indi- 
cates that decreases in MNase protection in a tupl A 
are associated with stress response pathways whereas 
increases in MNase protection are instead associated 
with translation and ribosome biogenesis pathways 
(Figure 6D-F and Supplementary Figure SI 3). 
Consistent with this, regulatory factors associated with 
both stress response and protein biosynthesis pathways 
show distinct Tupl -dependent chromatin structures over- 
lying their binding sites in vivo (Figure 7). The effects of 
Tupl regulation on chromatin structure at protein bio- 
synthesis genes appears to be indirect (Supplementary 
Figure SI 4), however, multiple lines of evidence are pre- 
sented which indicate that Tupl regulates a distinct 
subset of these genes, as the characteristic Tupl 
promoter structure is present for only the Tupl -bound 
fraction of protein biosynthesis genes (Supplementary 
Figures S13 and S14). 

Tupl has also been shown to regulate TF-binding spe- 
cificity previously, whereby a nucleosomal hair-trigger 
blocks the conditional binding of a TF until the appropri- 
ate environmental signal is presented (41). Therefore, we 
attempted to identify the most likely TFs that could 
respond to Tupl -regulation through chromatin by ex- 
amining the change in MNase protection centered on 
TFBSs (Figure 7). While our approach appears promising 
and identified the stress induced factor Xbpl and the 
general stress response factors Msn2 and Msn4, this 
analysis is limited because the majority of the data for 
the bound TFBS were obtained in experiments preformed 
during growth in rich media (50). Further experiments 
with the appropriate environment for conditional activity 
of suspected TFs will be required to identify conditional 
TFBSs whose binding may be regulated by Tupl- 
dependent nucleosome occupancy. Additionally, Tupl 
regulation of binding through chromatin will need to be 
confirmed by alternate approaches. However, from this 
analysis, it is clear that Tupl -dependent nucleosome oc- 
cupancy specifically overlies the binding sites of many 
stress response factors and its own recruiters, presenting 
another link between chromatin structure regulation by 
Tupl and stress-response pathways. 

Chromatin remodeling by Tupl represses transcription in 
conjunction with Isw2 

Tupl-Ssn6-mediated repression is believed to operate 
through a multistep mechanism involving HDAC re- 
cruitment, interactions with hypoacetylated histones and 
transcription machinery, and cooperation with other chro- 
matin remodeling complexes such as Isw2 (26). Combined 
mutations for genes encoding each of these factors have 
demonstrated larger effects on repression at individual loci 
in vivo, compared to single mutations, suggesting that co- 
operation between repressive factors enables separate, 



synergistic contributions to the regulation of genes (31). 
Moreover, unequal sensitivities of Tupl -regulated genes 
to the inactivation of different repression pathways 
in vivo also indicates that responses to Tupl repression 
may vary in a gene-specific fashion, whereby different 
groups of genes have evolved different strategies for 
utilizing Tupl-Ssn6 repression mechanisms (28). Here, 
we show that a majority (52%) of Tupl -remodeled pro- 
moters are also Isw2-remodeled (Figure 8A), indicating 
that the majority of Tupl -remodeling may be caused by 
its interaction with Isw2. Results extend individual loci 
studies defining a Tupl gene repression mechanism 
through Isw2-dependent nucleosome positioning (33,60) 
and demonstrate that this mechanism is functioning 
at many Tupl -regulated sites across the genome. 
Interestingly, an analysis of the association of site-specific 
Tupl recruiting cofactors binding with Tupl -remodeling 
events demonstrated that nearly all Tupl recruiting cofac- 
tors are associated with Tupl's activity on chromatin and 
that no single factor is solely responsible for directing 
Tupl's chromatin remodeling efforts. Additionally, the 
binding of Tupl, which follows the aforementioned re- 
cruiting efforts, is also not significantly predictive of 
which regions of chromatin Tupl selects to remodel. 
Therefore, it is clear from our analysis that the specificity 
of Tupl chromatin remodeling is not defined by Tupl 
recruitment or binding alone and requires additional 
cellular inputs. Further examination, focusing on 
post-translational modifications of Tupl -associated cofac- 
tors and their corresponding signaling pathways should 
shed light on how certain sites in the genome can have 
their chromatin remodeled by Tupl while other sites 
remain unaffected. 

Tupl targets open promoters to enhance downstream gene 
expression plasticity 

Despite the variable roles Tupl can play in regulating 
gene expression, Tupl -bound genes assume a distinct 
promoter packaging, consisting of low nucleosome occu- 
pancy/poorly positioned nucleosomes and a significantly 
wider NDR region (Figure 3C). Both the low nucleosome 
occupancy and wide NDR seen at Tupl -bound promoters 
in our study may reflect an abundance of easily digested 
nucleosomes at these loci that were digested away by our 
chromatin preparation or could reflect a more open chro- 
matin structure overall. Regardless of the possibilities, it is 
clear from our analysis that Tupl -bound genes are 
structured significantly differently than the rest of the 
yeast genome. This distinct promoter packaging also 
includes a wider distribution of conserved TF-binding 
sites and these loci are targeted by a larger number of 
/raws-activating factors in vivo. Together, this data dem- 
onstrates that Tupl regulation requires inputs from more 
regulatory factors at distinct promoter locations 
compared to other transcriptional regulatory mechanisms 
active in 5. cerevisiae. 

Overall, the characteristic chromatin signature of 
Tupl -regulated genes is consistent with previous studies 
of chromatin structure identifying distinct nucleosome oc- 
cupancy signatures in relation to transcriptional processes 
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(49), but what does this mean? Why does Tupl selectively 
target open regulatory regions? From analysis in Figure 3C, 
we can see that Tupl is not responsible for creating the 
distinct architecture at the genes it targets because deletion 
of Tupl did not revert the Tupl -bound promoters to a 
narrower NDR promoter with increased nucleosome oc- 
cupancy and improved positioning. Deletion of Tupl 
instead caused the opposite effect, producing a significant 
decrease in MNase protection at the —1 and —2 promoter 
nucleosomes. Therefore, knowing that this distinctive 
chromatin structure is not reliant on Tupl activity, we 
instead propose that Tupl is targeted based on its 
unique ability to enhance downstream expression plasti- 
city of these distinctly structured genes (Figure 5). We 
propose that the ability of Tupl to bind to both active 
and repressed genes (26), the diverse DNA-binding 
partners solicited by Tupl (42), the different regions of 
the Tupl-Ssn6 complex capable of interacting with 
DNA and protein partners (70), as well as the multiple 
mechanisms of transcriptional repression that Tupl can 
elicit, all enable Tupl to take advantage of these open pro- 
moters to promote downstream transcriptional plasticity. 
From our analysis, it seems likely that Tupl -dependent 
nucleosome occupancy is a key part of these transcription- 
al responses, as distinct changes in chromatin structure 
are seen at Tupl -regulated promoters depending on the 
associated gene expression program (Figure 6A-C), bio- 
logical process (Figure 6D-F), or underlying regulatory 
factor-binding sites (Figure 7) involved. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online. 
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